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(Dated: June 1, 2010) 

We derive the expressions for the anisotropic heterogeneous tensors of permittivity and perme- 
abihty associated with two-dimensional and three-dimensional carpets of an arbitrary shape. In the 
former case, we map a segment onto smooth curves whereas in the latter case we map a non convex 
region of the plane onto smooth surfaces. Importantly, these carpets display no singularity of the 
permeability and permeability tensor components, and this may lead to some broadband cloaking. 

PACS numbers: 42.70.Qs, 78.20Ci 



In 2006, the physicists Pendry, Schurig and Smith the- 
orized that a finite size object surrounded by a spherical 
coating consisting of a metamaterial might become in- 
visible for electromagnetic waves [1]. This is somewhat 
analogous to the alternative route to invisibility using 
conformal mappings (in the complex plane), preferred 
by Leonhardt [2] . These two proposals have captured the 
imagination of scientists working in the area of metama- 
terials. However, the former is not restricted by small 
wavelengths, and it has been experimentally validated in 
the microwave regime using a two-dimensional setup [3] . 

The underlying idea behind the cloaking using trans- 
formation optics is to map a point in optical space onto a 
spherical (invisibility) region. Back in 1984, the mathe- 
maticians Kohn and Vogelius noticed that one could find 
the conductivity of an object from static measurements 
on its boundary [4]. In the same vein, Greenleaf, Las- 
sas and Uhlmann looked in 2003 at an inverse problem 
where the Dirichlet to Neumann map defining a coating 
had the required properties to make a small conducting 
body nearly invisible |5]. But the important physical 
consequences had not be drawn by the mathematicians. 

Many authors have since then dedicated a fast grow- 
ing amount of work to the invisibility cloaking problem. 
Interestingly, there are alternative approaches, including 
some which make use of plasmonic properties of coated 
cylinders O |7j . These latter proposals are sometimes re- 
ferred to as external cloaking. The main advantage over 
the transformation optics approach is that there is no 
requirement for anisotropic heterogeneous permittivity 
and permeability, which is a consequence of the change 
of coordinates [5HTT]. However, external cloaking is nar- 
rowband in nature, whereas transformation optics allows 
for broadband cloaking, and works even in the intense 
near field limit when a source is located a couple of wave- 
lengths away from the cloak [121. Transformation optics 
can also be used to design generalized perfect lenses [13] . 

A severe limitation in the design of invisibility cloaks 
via transformation optics is the singular behaviour of the 
material parameters at the cloaks' inner boundary, which 
is a consequence of tearing apart the metric when one 
makes a hole in optical space (known in mathematics as 



blow up theory [5]). Physically, light has to curve its tra- 
jectory around the hole (or 'invisibility region'); Hence, 
to match the phase of a wave propagating in homoge- 
neous space, it must travel faster. One way to avoid 
such paradoxes is to approximate the cloaks' parameters 
using a homogenization approach, which leads to nearly 
ideal cloaking [TMTS] . Attractive theoretical proposals to 
avoid the cloaks' singularities include the design of nearly 
ideal (non-singular) two dimensional cloaks from a pro- 
jection of three dimensional ideal (but singular) cloaks 
[T7HI9] . An alternative route is to use a one-to-one map- 
ping to design an invisibility carpet, which is the bottom 
line of the bold proposal by Li and Pendry to conceal 
an object that is placed under a curved reflecting sur- 
face by imitating the reflection of a flat surface [2^ . The 
present letter is the first report of arbitrarily shaped two- 
dimensional and three-dimensional carpets. 

In electromagnetism, a change of coordinates induced 
by a geometric transform leads to the design of complex 
materials. For instance, if we start from a homogeneous 
and isotropic dielectric medium described by a permittiv- 
ity £ and a permeability jJ, = 1 (no magnetism), we end 
up with an inhomogeneous anisotropic material described 
by a transformation matrix T (also known as metric ten- 
sor) [Si [9j [m [12] . The permittivity and permeability in 
the transformed coordinates are now given by: 



eT" 



and ii' = /^T"^ where T = J^J/dei(J) 



(1) 
where J is the Jacobian matrix of the transformation. 
Importantly, we note that this material is magnetic. 

We now want to apply this recipe to design two- 
dimensional and three-dimensional carpets. Let us first 
consider the linear geometric transform: 

x' — X , a < X < b , 

"' ^'-^^^^^y + yi{^) , 0<y <y,ix) , (2) 
= z , —00 < z < -l-oo , 

where y' is a stretched vertical coordinate. It is easily 
seen that this linear geometric transform maps the seg- 
ment (a, b) of the horizontal axis y = onto the curve 



y' = Vi{x), and it leaves the curve y = 2/2 (a^) unchanged. 
Importantly, there is a one-to-one correspondence be- 
tween the segment and yi. The curves j/i and 1/2 are 
assumed to be differentiable, and this ensures that the 
carpet won't display any singularity on its inner bound- 
ary, as we shall now derive. 

The linear transform (pi) is expressed in a Cartesian 

/ 1 0~ 
basis as: S^x' = s '^ I where a = (1/2 - yi)lyi 

\ 1 
and from the chain rule 



dy_ 

dx' 



y2 



y' - V2 dyi 



{y2-yi)' 



dx 



2/1 



y' " yi dy2 



(2/2 - yiY 



dx 



(3) 



This leads to the inverse symmetric tensor T^^ which 
is fully described by five non vanishing entries in a Carte- 
sian basis: 



1 



(r-i)n--,(r-i)i2 = (r-i)2i = 



dy 



OL OT 



(4) 



It is interesting to look at the behaviour of the eigenval- 
ues of T~^ as these are the relevant quantities to compute 
the tensor components along the main optical axes: 



Al — — ; Ai 

a 



2a 



1 + a^ 



dy 
dx' 



+ (-irS/-4a2+(l + a2 + (0 



(5) 



We note that Ai and A^, i ~ 2,3, are strictly 



positive functions as obviously 1 



m 



'-4a2 



1 



and also a > 0. This 



establishes that T^^ is not a singular matrix for a two- 
dimensional carpet, which is a big advantage over two- 
dimensional cloaks obtained by blowing up a point onto 
a disc [U [12] : the transformation matrix is then sin- 
gular at the cloak's inner boundary (one eigenvalue goes 
to infinity, while the other two go to zero). 

For the sake of illustration, let us now consider a two- 
dimensional carpet that has inner and outer boundaries 
given by 



yi{x) = hi 



J 2\a I 



^ 7^ I +Ci sin di.h 



J 2 V c -' 



1 



(6) 
j = 1, 2, with hi = 0.2, /i2 = 0.4, ci = C2 = 0.01, di = 60, 
^2 = 50 and a = 0.3. 

We plot the profile of the three eigenvalues A^ along 
the inner and outer boundaries yi and j/2 of the carpet, 
as well as along the curve located half way from these, 
i.e. (2/2(2;) -I- yi{x))/2. We can see in Fig. fllthat none 



of the eigenvalues vanish and they satisfy the inequality 
< A2 < Ai < A3, a fact which can be also readily shown. 
We further note that Ai and A2 take values strictly within 
1.5 and 3.5. 

Thanks to the cylindrical geometry, the problem splits 
into p and s polarizations. In p polarization, we have: 



(7) 



V • e^ ViJa + Mo£ow>33^3 = 



in the carpet, where H; = {H3{x,y) — Hl{x,y))ez is 
the diffracted field parallel to the cylinder axis. Impor- 
tantly, H; satisfies the usual outgoing wave conditions 
as well as the Neumann data dH^/dn — dHl/dn on the 
ground plane and the inner boundary of the cloak, with 
H\{x,y)ez the incident field which is a beam generated 
by a constant field on a segment located at the upper 
left corner of the computational domain, and making an 
angle of 45 degrees with the horizontal axis. Moreover, 
e'rp is the upper left block diagonal part of e^ and /X33 the 

third diagonal entry of /u', as deduced from (fll) and (|4|. 

Such an anisotropic permittivity e^ could be achieved 

e.g. using some thin wires of metal diluted in dielectrics 
pn [55] to meet the condition that its eigenvalue A2 is 
lower than 1, see Fig. IT] Moreover, /igg — Ai involves 
some artificial magnetism which would require some res- 
onant elements such as split ring resonators [23j used in 
the design of the first invisibility cloak [3] . 



Hotof ?..=fl>JJ=(12,3) 




FIG. 1: Profile of the eigenvalues \i{x), i = 1, 2, 3 of T"^ on 
the inner boundary yi(x), the outer boundary 2/2(2;) and the 
centerline (5/2(3;) -I- yi{x))/2 of the carpet. 

In what follows, we consider a plane wave and a beam 
incident from the top at the wavelength A = 0.06. In Fig. 
[2] we report some computations where the plane wave is 
coming from above and the beam is incident from the 
top left corner, making an angle of = 45 degrees with 
the normal to the ground plane. We note that the field 
diffracted by the flat ground plane with infinite conduct- 



ing condition i.e. a mirror, cf. Fig. [2[d), and by an in- 
finite conducting object i.e. a curved mirror surrounded 
by the carpet, cf. Fig. l2jc), indeed superimpose. Of 
course, tlie field diffracted by the curved mirror on its 
own, cf. Fig. [2Fb), is much different. 

We then repeat the same simulation with a Gaussian 
beam in order to further analyse the effect of the carpet 
in a more realistic physical situation. We report these 
computations in Fig. |3] where it should be noted that 
the beam reflected by the carpet appears to have a waist 
closer to that of the incident beam than in the case of a 
flat mirror. This might be attributed to the fact that the 
optical path followed by the center of the beam is smaller 
in the case of a carpet. 
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FIG. 2; Diffraction by a plane wave and a beam at wavelength 
A = 0.06: we set H\ = 1 on the upper left side of the inner 
trapezoidal domain; 2D plot of the real part of the component 
Hz of the magnetic field, (a) Deformed mirror with a carpet 
under normal incidence; (b) Deformed mirror under oblique 
incidence; (c) Same as (b) with a carpet; (d) Flat mirror under 
oblique incidence. 

Let us finally consider the linear geometric transform: 

x' = x{r, 6*) , < r < p{B) . < 61 < 27r , 
y' ^y{r,9),O<r<p{0), < 9 < 27t , 

^' =^^^^'^£^^z + zi{x,y),0<z<Z2{x,y), 

(8) 
where z' is a stretched vertical coordinate. It is easily 
seen that this linear geometric transform maps the arbi- 
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FIG. 3: Diffraction by a Gaussian beam under oblique 
incidence at wavelength A — 0.06: we set Hg = 
ex.p{~l/2{l/V2{x + y)f /O.l^) on the left side of the inner 
triangle; 2D plot of the real part of the component H3 of 
the magnetic field, (a) Flat mirror; (b) Deformed mirror; (c) 
Same as (b) with a carpet. 



trary domain D = \J{u,v)e{oSf{i^('^''")'yi'^^'"))} within 
the plane xy onto the surface z' = zi{x,y), and leaves 
the surface z — Z2{x,y) unchanged. Importantly, there 
is a one-to-one correspondence between the domain D 
and the surfaces z' = Z\ and z' = zi- The surfaces z\ 
and zi are assumed to be differentiable, and this ensures 
that the carpet won't display any singularity on its inner 
boundary. 

The linear transform Q is expressed in a Cartesian 

/I 0' 

basis as: ^xx' = \ ^ 10 
\ dz dz 1 
\ dx' dy' OL 

and from the chain rule 



where a= {z-i — zi)/z2 



dz 
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dz 
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?■ dx 
dz\ 
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z\- 



z' — Z\ dz2 

iz2~ Zi 
z — Z\ 



?■ dx 
dz2 



(9) 



{z2-z^ydy {z2-ziydy 



This leads to the inverse symmetric tensor T ^ which 
is fully described by seven non vanishing entries in a 



Cartesian basis: 



ox' 




(10) 
Wc note that the entries of the transformation matrix 



dz 



in ( 10 ) reduce to those of \A\ when #^ vanishes. The 



corresponding eigenvalues have the similar structure to 
([5]) and are once again strictly positive and bounded, see 
Fig. |4j hence the material parameters are non-singular. 
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FIG. 4: (a) 3D plot of Ai; (b) Corresponding iso-contours; 
(c-f) Iso-contours of the eigenvalues \i{x,y), i = 2, 3 of T^^ 
given by (10 1 on the inner boundary Zi{x,y) (b,c,d) and the 
outer boundary Z2{x,y) (b,e,f) of the carpet; We note that 
iso-contours of Ai are the same whatever the altitude z. 

Let us now compute the total electromagnetic field for 
a plane wave normally incident upon a three-dimensional 
carpet. We have implemented the weak form of this 
scattering problem in the finite element package COM- 
SOL using second order finite edge elements which be- 
have nicely under geometric changes. Perfectly Matched 
Layers (PMLs), which can be seen as a stretch of coordi- 
nates, further enable us to model the unbounded domain. 
We choose the electric field E = {Ei,E2,E3){x,y, z) as 
the unknown in the Hilbert space //(curl, J7) = {v e 



[L^{n)] ,curlv e [L'^in)] } of curl-conforming fields [9], 
and therefore look for solutions of 



Vx 



/i' V X E 



- fc^e^E = 



(11) 



where fco = u^^qEq — lo/c is the wavenumber, c being 
the speed of light in vacuum, and e^ and //' are defined by 

Eqs. ([1]). Also, E = E^-I-Ed, where E^ is the incident field 
(here a field approximating a plane wave incident from 
the top which is generated by a constant field on a flat 
surface on the upper part of the computational domain) 
and Ed is the diffracted field which decreases inside the 
PMLs. We note that we also used this setting to retrieve 
our former computations assuming an electric field with 
the form {Ei^ _E2, 0) in ( 11 1 to take advantage of pull-back 
properties of edge-elements, leading again to Fig. |2]when 
we compute the curl of the numerical solution and plot 
the real part of the longitudinal component of (0, 0, H^) 
In this three-dimensional setting, we consider a plane 
wave incident from above at normal incidence: E^ = 
g-ifczg^^ with wavenumber k = 27r/0.3. The carpet has 
inner and outer surfaces given by: 



Zi{x,y) = hi [e 2(''^ ) - ij 
+c, sin ( di.h. 



e -2^ a ) 



1,2, 



(12) 



where p{0) = r(l — 0.1cos(56')) with r 



-- 0.2, h2 



2 and 
0.3, 



B — 2 arctan(?;/(a; + \J x"^ +2/^)) with h\ 

ci = 0.003, ci = 0.005 rfi = ^2 = 200 and a = 0.3. 

It is clearly seen from panels (c) and (d) in Fig. [5] that 
although of a complex non-convex shape, see (e), the car- 
pet refiects the plane wave nearly like a flat ground plane 
would. When the bump is not covered by the carpet, the 
scattering is much worse, see (a) and (b). 

Finally, we repeat these simulations for a Gaussian 
beam in oblique incidence (making an angle 7r/4 with 
the vertical axis) . This requires a computational domain 
shaped as a prism, see Fig. [6J We note that the plots of 
the field are indeed symmetric with the xOz plane in the 
case of a fiat mirror, see Fig. p[a-b), and a deformed mir- 
ror surrounded by the carpet, see Fig. ro[d-e). However, 
the diffraction by a deformed mirror is clearly giving rise 
to an asymmetric field, see Fig. l6lc). 

In this paper, we have shown that it is possible to 
design two-dimensional and three-dimensional carpets of 
an arbitrary shape, using a similar approach to Fourier- 
based cloaks [24] . Such carpets do not exhibit any singu- 
lar material parameters on their inner boundary, unlike 
invisibility cloaks, as they are based upon a one-to-one 
geometric transform. The next step towards the realiza- 
tion of such carpets might involve some structural ele- 
ments such as conducting thin-straight wires and split 
ring resonators [13] to tune the permittivity and perme- 
ability to required values depending upon light polariza- 
tion. The rapid experimental progress in the construc- 
tion of carpets getting close to optical frequencies 
suggests that our designs might soon come to life 
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FIG. 5: DifTraction of a normally incident plane wave by a 
deformed mirror surrounded by a 3D carpet at wavelength 
A = 0.3; (a)-(d) 2D plots of the real part of the component 
E3 of the electric field in the planes xOz, yOz for the bump on 
its own (upper panel) and with a carpet (middle panel); (e) 
3D plot of the real part of E3: Cartesian Perfectly Matched 
Layers were implemented in the domains surrounding the cen- 
tral cubic region (lower panel). 
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